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The MACHO Project is searching for galactic dark matter in the form of massive compact halo objects (Machos) . 
Millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge are 
photometrically monitored in an attempt to detect rare gravitational microlensing events caused by otherwise 
invisible Machos. Analysis of two years of photometry on 8.5 million stars in the LMC reveals 8 candidate 
microlensing events, far more than the ~ 1 event expected from lensing by low-mass stars in known galactic 
populations. From these eight events we estimate the optical depth towards the LMC from events with 2 < t < 200 
days to be r| 00 w 2.9^J'g x 10 -7 . This exceeds the optical depth of 0.5 x 10 -7 expected from known stars and is 
to be compared with an optical depth of 4.7 x 10~ 7 predicted for a "standard" halo composed entirely of Machos. 
The total mass in this lensing population is w 2^J'y x 10 11 Mq (within 50 kpc from the Galactic center). Event 
timescales yield a most probable Macho mass of 0.5^0 2 Mg, although this value is quite model dependent. 



1. INTRODUCTION 

Galactic dark matter in the form of Machos can 
be detected by its occasional gravitational mi- 
crolensing of extragalactic stars [ Q g, g|. This 
effect occurs when a Macho of mass m is in 
close alignment with a background star and acts 
as a gravitational lens, magnifying and distort- 



ing the stellar image. Microlensing refers to the 
situation in which image distortion is not de- 
tectable (<C 1 arcsec) and the only visible ef- 
fect is an apparent amplification, A = (u 2 + 
2)/u\/u 2 + 4. Here u — bjrE is the impact pa- 
rameter in units of the Einstein radius te = 
y/ AGttiDolDls/ c 2 D s where D l, D ls & D s 
are distances between observer, source and lens[ 
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|40. This effect is transient with a timescaie, i = 
2?\e /v± , set by the mass and iocation of the lens 
and its motion relative to the line of sight. For a 
source in the Magellanic Clouds and a lens at 10 
kpc, r E w &y?m/MQ AU and i w 14O^/m/M 
days for v± = 200 km s _1 . With ^daily observa- 
tions over periods of years, Machos in the range 
10~ 4 — 10 2 M0 are experimentally accessible. 

Because the "area" of the lens is proportional 
to its mass, the fraction of solid angle occupied 
by lenses, r, depends only on the mass density 
of Machos and not on their particular sizes. Ap- 
plication of the virial theorem shows that t ~ 
(v c /c) 2 ~5 x 10 -7 for a Macho-dominated halo. 
Here v c is the circular velocity about the cen- 
ter of the galaxy. Thus detection of these very 
rare events requires many repeated observations 
of millions of sources. The optical depth towards 
the Magellanic Clouds due to a Macho-dominated 
halo still comfortably exceeds that from known 
populations of low-mass stars, so that microlens- 
ing surveys directly probe the Macho content of 
the halo. 

Several groups are now engaged in large scale 
photometric surveys to search for dark matter 
in the form Machos using the high density of 
stars in in the Magellanic Clouds, M31 and the 
Galactic bulge as a backdrop for gravitational 
microlensing. Four such surveys, DUO, EROS, 
MACHO and OGLE, have detected in excess 
of 100 events that appear to be microlensing [ 

1 §, || 6 [I of which ~ 10 are to- 
ward the LMC, the remainder being toward the 
Galactic bulge. Most of these events have charac- 
teristic timescales t ~ 20 — 60 days. Searches for 
short-timescale events with 1 hour < t < 10 days 
have revealed no candidates to date and set in- 
teresting limits on low-mass Machos [ [l 
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Detailed reviews of microlensing are given by [ |l7fl 
and [plij|; a basic outline of the MACHO project, 
and analysis procedures is given in [[H (hereafter 
A96) and a more detailed analysis of this sample 
is found in [||. 

2. OBSERVATIONS & PHOTOMETRY 

The MACHO Project has had full-time use of 
the 1.27-meter telescope at Mount Stromlo Ob- 
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Figure 1. An R-band image of the LMC, 8.2 de- 
grees on a side (G. Bothun, private communica- 
tion), showing the locations of the 22 MACHO 
fields used here. Event locations are are shown 
with circles. 



servatory, Australia, since mid 1992 in an ex- 
tended run until the end of the century. Details 
of the telescope system are given in [ 20 and of 
the camera system in [ |2l], ^2j . Briefly, corrective 



optics and a dichroic beam splitter are used to 
simultaneously image a 42 x 42 arcmin 2 field in 
two colors, using two mosaics of four 2048 2 pixel 
CCDs. 

As of this writing, over 44000 exposures have 
been taken with the system, over 3 TBytes of 
raw image data. About 60% are of the LMC, 
the remainder are towards the Galactic bulge and 
SMC. In this paper, we consider only the first 2.1 
years of data from 22 well-sampled LMC fields, 
located in the central ~ 5° x 3° of the LMC and 
shown in Figure [IJ 

The observations described here comprise 
10827 images distributed over the 22 fields. These 
include observations over a time span of 840 days 
from mid 1992 to late 1994. The mean number 
of exposures per field is 492, with a range from 
300 to 785. The sampling varies between fields 
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Table 1 

Event selection criteria 





Description 


Value 


1 


Coverage 


t < 300 days & t max £ data range 


2 


Constant baseline 


xli-ouJNdof < 4 with N dof > 40 


3 


Amplification 


A max > 1 + 2(7 


4 


Significance 


6 measurements > ltr high 


5 


Variable (Bumpers & very red stars) 


V>17.5 & V-R<0.9 


6 


SN 1987A light echos 


10'xlO' region excluded 


7 


Peak fit 


AxV^cakMiof) > 200 


8 


Crowding 


> 95% of measurements have low crowding 


9 


Crowding 


/ crd < 701og 10 [Ar7(xL/AM] - 45 


10 


Significance 


Ax7(xL/^dof) > 500 


11 


Amplification 


^4max > 1.75 



since the higher priority fields were often observed 
twice per night with a ~ 4 hour spacing for sen- 
sitivity to very short timescale events [ [l5|, M . 

The photometric reduction procedure here is 
identical to that described in A96; briefly a good- 
quality image of each field is chosen as a template, 
and used to generate a list of stellar positions and 
magnitudes. All other images are aligned with 
the template using fiducial stars, and a PSF is 
determined from these. Then, the flux of all other 
stars is fitted using the known positions and PSF. 
For each measurement we also compute an error 
estimate and six quality flags; these flags are the 
object type, the x 2 of the PSF fit, a crowding 
parameter, a local sky estimate, and the fraction 
of the star's flux rejected due to bad pixels and 
cosmic rays. Suspect measurements are removed 
based on these quality flags and the resulting data 
are organized into stellar lightcurves. 

3. EVENT DETECTION 

The MACHO Project is the largest survey of 
astronomical variability in history and as such 
encompasses many previously unknown back- 
grounds to microlensing. The data used here 
comprise some 9 billion photometric measure- 
ments. Although we use objective selection cri- 
teria (cuts) to discriminate genuine microlensing 
from stellar variability and systematic photome- 
try errors, these cuts are developed after extensive 
examination of the MACHO database. However, 



as long as the experiment's event detection effi- 
ciency is calculated properly and the selection 
criteria are sufficiently stringent to accept only 
real microlensing events, changes in the selection 
criteria will be accounted for in the efficiency cal- 
culations, and the details will not bias the final 
results. 

Our experience with the Galactic bulge has 
been helpful in developing the current cuts. Since 
our previous analysis of a smaller set of LMC 
data (A96), we have analyzed a large set of data 
towards the Galactic bulge [ [| |l(j, which has 
yielded over 80 microlensing events. These obser- 
vations provide a more realistic set of microlens- 
ing events than the artificial events we previously 
used to test our analysis procedures. A number 
of these events fail the original cuts because they 
are not well described by the "normal" microlens- 
ing light curve which assumes a single point lens, 
constant velocities, and an unblended source star 
but nevertheless appear by eye to be very high 
quality microlensing events. In contrast events 2 
and 3 of A96 appear far less striking. 

This has motivated us to modify our selection 
criteria, shown in Table |l|, to have less strin- 
gent requirements on the functional form of the 
lightcurve. We compensate for this by demand- 
ing higher signal to noise i. e. higher amplification 
and greater statistical significance. 

Event detection proceeds in several stages. In 
the first stage, the lightcurves of §^ are con- 
volved with 7, 15 and 45 day matched filters as in 
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A96 and low level triggers are fit with the "stan- 
dard" microlensing curve A(u(t)), where u(t) = 

"J w min + ~ tmax) 2 /i 2 - We then compute for 
each light curve a set of statistics and apply the 
cuts shown in Table |l| to distinguish microlensing 
from the background. A measure of significance 
is A x 2 /{x 2 m i/N doi ), where A x 2 = xlonst ~ X 2 mi is 
the improvement in \ 2 between a constant bright- 
ness fit and a microlensing fit. 
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Figure 2. The final cuts for selection of microlens- 
ing candidates. The y-axis is the fit A max . The 
solid lines show the final cuts 10 and 11; the cir- 
cles in the upper right are the 12 events (10 stars) 
passing all cuts. 



Two of the principal cuts are summarized in 
Figure ||, which illustrates cuts 10 and 11 for all 
events that pass cuts 1-4 & 6. Events that pass 
all cuts are indicated with circles, while the other 
open symbols and crosses indicate events which 
failed 1 or more of cuts 5 & 7-9.. Open squares 
indicate events which fail cut 5 on magnitude, 
and open triangles indicate events which fail cut 
7 on peak fit. The crosses indicate events which 
fail cuts 8 or 9. The lightcurves passing all of 
these selection criteria are further investigated, 



as outlined below. 
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Figure 3. Candidate lightcurves 1 & 4 



4. MICROLENSING CANDIDATES 

Twelve lightcurves pass the cuts discussed 
above, indicated by the large open and filled cir- 
cles in the upper right hand region of Figure |^. 
Some lightcurves are shown in Figures [| & || and 
event fit parameters are listed in Table |^. Four 
of these lightcurves (la, lb, 12a and 12b) cor- 
respond to stars which occur in field overlap re- 
gions; the two lightcurves for each star are based 
on independent data and reductions and are de- 
tected independently. 

Event 1 was our first discovery [ and has 
remained constant in following years. The star's 
spectrum is that of a normal clump giant at the 
radial velocity of the LMC [ ^3) . 

The low signal-to-noise candidates 2 and 3 from 
A96 do not pass the final cuts used for this data 
set and are not included in following calculations, 
but we have numbered the present set of 10 candi- 
dates 1, 4 ... 12 to avoid possible ambiguity. The 
EROS group has informed us that the star in- 
volved in event 2 may have brightened in 1990, 
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Figure 4. Candidate lightcurves 5-7 



background galaxy and is almost certainly a su- 
pernova in that galaxy. Event 12 may in fact be 
microlensing. Our efficiency analysis, described 
below, does not take into account this "ampli- 
fication bias" therefore it is necessary to discard 
both of these events for the following calculations. 

These events, as shown in Figure [j], do not clus- 
ter spatially or in any particular stellar popula- 
tion. In addition, the distribution of fit A max is 
consistent with the expected uniform distribution 
of Wmin- We classify events 1,4,5 & 9 as 'excellent' 
microlensing candidates, events 6-8 as 'good' and 
event 10 as 'marginal'. Events 11 & 12 must be re- 
jected as explained above, and events 2 & 3 from 
A96 do not pass the current cuts. We define two 
sets of events for the following analysis. The first 
set contains all eight events passing the cuts with- 
out amplification bias. The second, more conser- 
vative, set of 6 events eliminates event 9 as a pos- 
sible LMC lens and event 10 as a possible variable 
star. 



and it has also shown indication of further bright- 
ening episodes in our subsequent data; thus it is 
probably a variable star. 

Event 4 was the first LMC microlensing candi- 
date detected in progress by our real-time Alert 
system [ |l0| and announced in IAU Circular 6095. 
This event is confirmed by additional photome- 
try and spectra, obtained near peak amplifica- 
tion! |§. 

Event 9 passes our objective data cuts, but it 
appears to be lensing by a binary system. This 
event may resolve the extended size of the source 
star, providing a measure of the relative proper 
motion of the lens and source, which indicates the 
lens may reside in the LMC [ ||[| . 

Event 10 passes all of our cuts, but our experi- 
ence suggests that it may be a variable star. The 
asymmetry seen in the light curve, a rapid rise 
with with a more gradual fall, is typical behavior 
one expects from variable stars. 

Events 11 and 12a/b are indicated with open 
circles in Figure || as they are only present in 
our reference template images because of their 
amplification during that particular observation. 
They are no longer detectable in images that go 
somewhat deeper. Event 11 is superimposed on a 



5. DETECTION EFFICIENCY 

The detection probability for individual events 
depends on many factors, e.g. the 3 event pa- 
rameters A max ,t,t max , and the unlensed stellar 
magnitude, as well as our observing strategy and 
weather conditions. However, we can average 
over the known distributions of A max and t max 
and stellar magnitude, using the known time- 
sampling and weather conditions, to derive our 
efficiency, shown in Figure as a function only 
of event timescale £ (t). 

We have computed our detection efficiency 
using an essentially identical method to that 
outlined in A96, simply generating simulated 
microlensing events with t logarithmically dis- 
tributed in the range 0.3-1000 days and adding 
these simulated events into the extended time- 
span of observations. The Monte-Carlo proce- 
dure takes into account the actual spacing and 
error bars of the observations, so any variations in 
sampling frequency, weather, seeing etc. between 
the first and second year data are automatically 
accounted for. 

Efficiency is defined as the fraction of input 
events with u m j n < 1 which pass our cuts; since 
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Table 2 



Event parameters & contribution to r 



Event 


V 


V-R 


^max 


t 




X 2 /N doi 


tbi 


n 


la 


19.6 


0.6 


57.08(3) 


34.7(3) 


7.2(1) 


1.420 


38.8 


1.8 x 10~ 8 


lb 


19.6 


0.6 


57.26(4) 


34.3(3) 


7.5(3) 


1.134 


4 


20.0 


0.2 


647.2(2) 


46(2) 


3.00(4) 


1.416 


52 


2.3 x 10~ 8 


5 


20.7 


0.4 


24.0(3) 


82(2) 


58(5) 


1.680 


88 


3.5 x 10~ 8 


6 


19.6 


0.3 


197.5(7) 


87(4) 


2.14(4) 


0.873 


100 


4.1 x 10~ 8 


7 


20.7 


0.4 


463.0(3) 


115(3) 


6.16(10) 


1.447 


131 


6.0 x 10~ 8 


8 


20.1 


0.3 


388.4(5) 


62(2) 


2.24(5) 


2.218 


70 


2.8 x 10" 8 


9 1 


19.3 


0.3 


603.04(2) 


143.4(2) 





1.755 


143 


6.6 x 10~ 8 


10 


19.4 


0.2 


205.3(3) 


42(1) 


2.36(5) 


1.982 


47 


2.1 x 10~ 8 


11 


21.5 


0.4 


-8.6(3) 


266(9) 


11.9(4) 


2.964 






12a 


21.2 


0.3 


-10.0(3) 


138(5) 


7.2(4) 


1.487 






12b 


21.2 


0.3 


-11.4(8) 


162(8) 


6.8(2) 


1.536 
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Figure 5. Microlensing detection efficiency (nor- 
malized to M m in < 1) for the 2-year Macho data, 
as a function of event timescale t. 



we use a cut of A maK > 1.75 or u min < 0.661, 
our efficiency is constrained to be less than 0.661. 
This efficiency is defined relative to an 'expo- 
sure' of E = 1.82 x 10 7 star-years, which arises 
from the 7.9 million well measured, independent 
lightcurves in our sample multiplied by the 840- 
day span of observations. 

We calculate our efficiency with two levels of 



realism. The first is the 'sampling' efficiency, in 
which we neglect stellar blending and assume that 
all the additional flux in a microlensing event is 
recovered by the photometry code. The second is 
the 'photometric' efficiency, where we add artifi- 
cial stars to a representative set of real images, re- 
run the photometry code to create look-up tables 
of added vs. recovered flux. These look-up tables 
are then used to generate artificial microlensing 
lightcurves in the same way as above. This 'pho- 
tometric' efficiency is more realistic, and is typi- 
cally ~ 20% lower than the 'sampling' efficiency 
for timescales less than 200 days. 

The photometric efficiency is based on all stars 
in our fields, even those not uniquely identified 
because of S/N or crowding effects. These are ac- 
counted for by integrating the detection efficiency 
per star over a corrected luminosity function (LF) 
as in A96. This corrected LF is truncated about 
one magnitude beyond where our measured LF 
becomes seriously incomplete. However, the real 
LF continues to rise as 10°- 5m for several magni- 
tudes beyond this cutoff so there should be an ad- 
ditional contribution to our exposure from these 
fainter stars. We have tried several different mag- 
nitude cutoffs and it appears that our exposure 
{EE) is converging near the cutoff used in this 
paper for events with t < 150 days. 

In the photometric efficiency Monte Carlo, we 
find that in addition to reducing our efficiency, 
extreme crowding tends to make events appear 
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lower amplitude and shorter. This effect is due 
to "blending" or the tendency of a group of sev- 
eral superimposed stars to be detected as a single 
object. We make statistical corrections for this 
effect in column 8 of Table ||. Our estimated sys- 
tematic error in the efficiency determination is 
~ 10%, less than our statistical error with eight 
events [ |||. 

The efficiency for timescales t ~ 10 — 60 days 
is lower than that from A96 by ~ 10% primarily 
because of the tighter A max and A% 2 cuts. For 
longer events, the two year data is more sensitive 
because of its longer time span. 

6. OPTICAL DEPTH 

The simplest measurable quantity in a gravita- 
tional microlensing experiment is the microlens- 
ing optical depth. In the the limit t«1, this is 
the instantaneous probability that a random star 
is magnified by a lensing object by more than 
a factor of 1.34 (it m i n < 1). This is related to 
the mass in microlensing objects along the line of 
sight to the source stars by 

c 2 J Dos 

Thus, it depends only on the density profile of 
lenses, not on their masses or velocities. Experi- 
mentally, one can obtain an unbiased estimate of 
the optical depth as 

T^ aS =--Y^J^-y (2) 

where E is the total exposure (in star- years), U 
is the Einstein ring diameter crossing time, and 
£(ii) is the detection efficiency. Here, and be- 
low, we use the blend corrected values, tu, from 
column 8 of Table § 

Confidence levels on T meas are determined from 
Monte Carlo experiments, in which the number of 
events is a Poisson variable with adjustable mean, 
and event timescales are drawn at random from 
the observed set. The mean is fixed so that 16% of 
experiments yield an optical depth larger/smaller 
than that measured. The mean optical depth for 
these experiments is then the "la" lower/upper 
bound on r mcas . 



Using our full sample of 8 events, we find an op- 
tical depth for events of duration 2 days < t < 200 
days of t| 00 = 2.9 +Jg x 10~ 7 . If we subtract the 
predicted background microlensing optical depth 
of Tbackgnd = 0.5 x 10~ 7 (from Table || below) , we 
find that the observed excess is about 50% of the 
predicted microlensing optical depth for a "stan- 
dard" all-Macho halo of e.g. [ [|J. Alternatively, 
we can estimate the optical depth due only to the 
halo by using the 6 event subsample defined in 
section § for which t| 00 = x 10~ 7 , about 

45% of the optical depth predicted by a "stan- 
dard" all-Macho halo. 

Uncertainties can be much larger than those 
calculated above if our measured distribution of 
timescales is a poor estimate of the real timescale 
distribution, especially if very long events are pos- 
sible. It is worth noting that errors are consider- 
ably larger than Poisson because they are dom- 
inated by the shot noise on the small number 
of long events rather than the total number of 
events. 

7. NON-HALO POPULATIONS 

An in depth assessment of the microlensing 
background is necessary to determine if r mea s 
and A mcas are significantly greater than that ex- 
pected from non-halo populations. As noted by 
[ |27], |2[|] and A96, low-mass stars in the 
Galaxy and LMC will give rise to some microlens- 
ing events. However, these authors find that the 
optical depth from known stars is only < 10% of 
that from an all-Macho halo. 

Shown in Table | are estimates of stellar lensing 
rates of A96, using the same model parameters 
for the thin and thick disks, the spheroid and the 
LMC disk, assuming a Scalo PDMF for all pop- 
ulations, and simply updating the total exposure 
for our 2-year dataset. In our present sample, we 
expect a total background of 1.14 events from all 
know stellar populations. 

A Poisson distribution with a mean of 1.14 
gives probabilities to observe A Q b s > 3,4,5,6,7 
as 10.7%, 2.9%, 0.6%, 0.11%, 0.02% respectively. 
Thus we see that if only three of our events are 
genuine microlensing, the evidence for an excess 
is modest, whereas if more than four of our events 
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Table 3 



Microlensing by stars 



Population 


r(10" 7 ) 


(t) (days) 


(Z)(kpc) 


r(io- 7 yr _1 ) 




Thin disk 


0.15 


112 


0.96 


0.62 


0.29 


Thick disk 


0.036 


105 


3.0 


0.16 


0.075 


Spheroid 


0.029 


95 


8.2 


0.14 


0.066 


LMC center 


0.53 


93 


49.8 


2.66 


(1.19) 


LMC average 


0.32 


93 


49.8 


1.60 


0.71 


Standard halo 


4.7 


89 


14.4 


24.3 


11.2 



are microlensing there is very strong evidence for 
an excess over stellar lensing alone. 

If one supposes that only the four "best" events 
result from lensing by stars and the other candi- 
dates are not microlensing (e.g. they are variable 
stars) then there is only marginal evidence for 
Machos in the halo. The flaw in this argument is 
that the distribution of microlensing magnifica- 
tions is given a priori, and the "best" microlens- 
ing candidates are preferentially the high magni- 
fication ones. We should expect to find a mixture 
of high and low magnification candidates, and it 
is very unlikely that all events would be high mag- 
nification given our detection efficiency. By this 
argument, the surplus of events is not critically 
dependent on the lower quality candidates and it 
is difficult to explain our results by stellar lensing 
alone. 

8. LIKELIHOOD ANALYSIS 

Since the timescale of a lensing event is pro- 
portional to \/m, we may use the observed event 
timescales to estimate the lens masses, given a 
particular halo model to fix the phase space distri- 
bution of lenses. We apply a maximum-likelihood 
method as in A96 using the "standard" halo 
where a fraction / of the dark halo is made of 
Machos with a unique mass m (the remaining 
1 — / of the halo is assumed to be unobserved dark 
matter). The resulting likelihood contours, as- 
suming a delta-function mass function, are shown 
in Figure [| the probabilities are computed us- 
ing a Bayesian method with a prior uniform in / 
and log m. The peak of the likelihood contours 
gives the most probable mass and halo fraction 
for this halo model. For the 6 event sample we 



find m 2 D = 0.41 M , and f 2D = 0.51. 

We calculate one-dimensional likelihoods by in- 
tegrating over the other parameter. For the six 
event sample this yields a most likely Macho mass 
m = 0.461°;?? M , and most likely halo fraction 
/ = 0-50±g;|g. The errors given are 68% CL. It is 
important to note the large extent of the contours 
in Figure l(| This is mostly due to the small num- 
ber of events. The 95% CL contour includes halo 
fractions from 17% to 100%, and Machos masses 
from 0.12 to 1.2 M . 

The most probable mass and halo fraction are 
both larger than our results from year-1, which 
were m = 0.06 M Q and / = 0.2, though there 
is a reasonable degree of overlap of the contours. 
The year-1 90% CL contour is shown as the light 
line in Figure ^. The most probable value of each 
analysis lies outside the 90% CL contour of the 
other analysis. This is primarily due to the up- 
ward shift in (t) . Events 2 & 3 from A96, dropped 
in this analysis, are both shorter than any events 
in the current set. The shift in (tj is also to be 
expected as our maximum efficiency shifts toward 
longer events with the increasing baseline of the 
experiment. 

9. DISCUSSION 

One natural explanation of the results pre- 
sented here is that a substantial fraction of the 
Galactic dark halo is in the form of Machos. The 
fact that the observed events have relatively long 
timescales suggests that (for canonical halo mod- 
els) the lenses have masses above ^0.1 M Q , with 
a most probable mass ~ 0.5 M . If so, they 
cannot be ordinary hydrogen-burning stars since 
there are strong direct limits on such objects from 
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Figure 6. Likelihood contours of Macho mass 
m and halo Macho fraction / for a delta- 
function mass distribution, for halo model S. 
The most likely value is indicated with a +, 
and the contours enclose total probabilities 
of 34%, 68%, 90%, 95%, 99%, using a Bayesian 
method as described in § || The light line shows 
the 90% contour from A96. 



counts of faint red stars (e.g. [ ^9|, |3(], [n|); thus 
stellar remnants such as white dwarfs appear to 
be an obvious possibility. 

There are some theoretical difficulties with the 
white-dwarf hypothesis (e.g. [ First, 
the initial mass function must be fairly sharply 
peaked between ^2 — 6 M Q to avoid overproduc- 
ing either low-mass stars (which survive to the 
present) or high-mass stars (which as type-II SNe 
overproduce metals). A second difficulty is that 
the high luminosity of the B and A progenitors 
may exceed the observed faint galaxy counts. 

However, white dwarfs are one of the few dark 
matter candidates which are known to exist in 
large numbers. Also, it has recently become clear 
(e.g. [ Q) that the mass of hot gas in rich galaxy 
clusters greatly exceeds that in stars, and further- 
more this gas is relatively metal-rich with an iron 
abundance ~ 0.3 solar. This might suggest that 



most of the baryons have been processed through 
massive stars, which have since died leaving a 
population of remnants and metal-rich gas. This 
scenario has been explored by [ |3^| , who suggest 
that it may be natural to have 40-100% of the 
Galactic dark matter in white dwarfs. 

The observational limits on the local density of 
white dwarfs are a strong function of their age [ 
p6[ [37J. Very recently, a limit from the Hubble 
Deep Field has been given by [ |3l| ; they find that 
white dwarfs with Mj < 16 contribute < 100% 
of the halo density, and those with Mi < 15 con- 
tribute < 33%. 

Our exploration of halo models also shows 
that for 'minimal' halos, the timescales may still 
be consistent with substellar Machos just below 
O.IMq. Also, models with a substantial degree 
of halo rotation may lead to smaller mass esti- 
mates, since the rotation could reduce the trans- 
verse component of the Macho velocities. Thus, 
brown dwarfs cannot be ruled out as yet. 

Finally, it is worth noting that microlcnsing 
is sensitive to any compact objects, irrespective 
of composition, as long as they are smaller than 
their Einstein radii. Thus, more exotic objects 
such as primordial black holes and strange stars 
are possible candidates that are virtually unde- 
tectable by direct searches. 

The most direct method for identifying the lens 
population is to measure microlensing 'parallax' 
by simultaneously observing events from earth 
and a small telescope in Solar orbit [ This 
measures the velocity of the lens projected to the 
Solar system, which provides a definitive proof 
of microlensing, and discrimination between disk, 
halo and LMC lenses on an event-by-event basis. 

To summarize, our results indicate a microlens- 
ing optical depth of « 3 x 10 -7 or a Macho mass 
within 50 kpc of « 2 x 10 11 M Q . This provides 
evidence that Machos with masses in the range 
0.05 — 1 M© contribute a substantial fraction of 
the galactic dark halo. Continued observations 
from this and other projects should clarify this in 
the next few years. 
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